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Abstract One of the key challenges in modeling the dy- 
namics of contagion phenomena is to understand how the 
structure of social interactions shapes the time course of 
a disease. Complex network theory has provided signif- 
icant advances in this context. However, awareness of 
an epidemic in a population typically yields behavioral 
changes that correspond to changes in the network struc- 
ture on which the disease evolves. This feedback mech- 
anism has not been investigated in depth. For example, 
one would intuitively expect susceptible individuals to 
avoid other infcctcds. However, doctors treating patients 
or parents tending sick children may also increase the 
amount of contact made with an infecteds, in an effort 
to speed up recovery but also exposing themselves to 
higher risks of infection. We study the role of these care- 
taker links in an adaptive network models where indi- 
viduals react to a disease by increasing or decreasing the 
amount of contact they make with infected individuals. 
We find that pure avoidance, with only few caretaker 
links, is the best strategy for curtailing an SIS disease in 
networks that possess a large topological variability. In 
more homogeneous networks, disease prevalence is de- 
creased for low concentrations of caretakers whereas a 
high prevalence emerges if caretaker concentration passes 
a well defined critical value. 



1 Introduction 



Physicists have taken numerous approaches to model- 
5-H ing infectious diseases, ranging from simple, determin- 
istic compartmental models that qualitatively describe 
disease dynamics in single populations [3] , to highly com- 
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plex, stochastic metapopulation models that can account 
for the spread of emergent infectious diseases on a global 
scale [HUTMTfj] . Simple models, designed to investigate 
the basic mechanisms underlying disease dynamics, typi- 
cally assume that a population is well-mixed, that inter- 
acting individuals are identical and that stochastic ef- 
fects are negligible [3l|9]. On the other hand, complex 
computational models are manufactured to predict the 
time-course of actual emergent infectious diseases such 
as H1N1 in 2009 [5], SARS in 2003 [19] quantitatively. 
They typically take into account data on social variabil- 
ity, age structure, spatial heterogeneity, seasonal varia- 
tion of disease dynamic parameters, multi-scale mobility 
networks, and account for stochastic effects. Both classes 
of models fulfill equally important, complementary, but 
almost mutually exclusive purposes. 

Theoretical epidemiology experienced a major thrust 
with the advent complex network theory and its intro- 
duction into the field [HI24] • The study of network prop- 
erties substantially advanced our understanding of dis- 
ease dynamic phenomena on multiple levels [26] . On one 
hand, networks were used as a model for inter-individual 
relationships (social networks) [53]. On the other hand, 
the network approach was applied on a larger scale, mod- 
eling mobility and transport between populations |191 

mi. 

The use of network theoretical concepts allowed re- 
searchers to investigate how topological properties of un- 
derlying networks shape the contagion processes that 
evolve on them [2^1[2Hl[71ll21l28j . In the context of epi- 
demiology, mapping structural features of networks to 
properties of the spread of the disease substantially in- 
creased the predictive power of models and our under- 
standing of epidemic phenomena. 

Although it is intuitive and plausible that network 
features determine the spread of a disease, it is equally 
plausible that an epidemic reshapes the structure of the 
underlying networks. For example, in response to in- 
formation on an ongoing epidemic, people may change 
their behavior. They may decide to wear face masks, 
avoid contacts, and travel less. Surprisingly, this feedback 
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mechanism has been neglected even in some of the most 
detailed and sophisticated modeling approaches [T41IT5]. 
Topological properties of social networks affect disease 
dynamics, and the disease then feeds back to change the 
topology of the network. In order to understand the dy- 
namics of contagion phenomena in a population, it is 
vital to understand the consequences of this feedback 
mechanism. 

Networks that change their structure in response to 
their environment are called adaptive [16,8,18,21 . In a 
recent study, Gross et al. proposed a simple adaptive 
network scheme, based on a rewiring rule, to understand 
how individuals' behavioral changes impact on the time 
course of an epidemic. In this model, susceptible individ- 
uals are allowed to protect themselves from infection by 
rewiring their existing links |17] . Specifically with prob- 
ability w a susceptible breaks the relationship with an 
infected person and forms a new link to another, ran- 
domly selected susceptible. Despite the simplicity of this 
approach, the mechanism can generate an abundance of 
interesting phenomena including hysteresis and multi- 
stability 

Although this mechanism is attractive, the response 
to an ongoing epidemic in a population has many facets. 
Not only do individuals avoid other infected individuals 
(negative response). In many scenarios, individuals in- 
crease their interaction with infected individuals (pos- 
itive response), particularly in hospital scenarios, and 
families in which individuals adopt the role of a care- 
taker. Potentially, these positive responses can facilitate 
disease proliferation in a population and yield a higher 
disease prevalence. However, caretaker activity can have 
a positive effect on infected individuals, for example by 
increasing a person's recovery rate. A key question is how 
these effects interact and under what circumstances care- 
taker activity has a net positive or negative effect and 
how these effects play out in different network topolo- 
gies. 

Here we propose and investigate these questions us- 
ing an adaptive network model. We consider two types 
of networks. First, the generic Erdos-Renyi random net- 
work with binomial degree distribution, where each pair 
of nodes is linked with constant probability per |13l 
l2~i] . We also consider Barabasi- Albert scale- free networks 
with power law degree distributions, which more closely 
mimic the heterogeneity in social interactions. Dynamics 
on scale-free networks have a number of important prop- 
erties. For instance, they lack epidemic thresholds and 
are immune to random immunization due to strong con- 
nectivity fluctuations |26l[25l [24 l l6 l l27] . Thus diseases on 
scale-free networks are difficult to avoid, and once they 
take hold, they are difficult to eradicate. We will show 
that in scale free topologies the highest disease extinction 
probabilities occur in the total absence of caretakers, a 
surprising result which suggests that caretaker relation- 
ships (including doctor/patient relationships) should be 
minimized in those systems. For Erdos-Renyi networks 



we observe a critical caretaker proportion which mini- 
mizes disease severity and beyond which additional care- 
takers increase disease prevalence. 



2 Model description 

We consider a network with a constant number of nodes 
N, representing individuals in a population. Each node is 
either susceptible (S) or infected (I). We denote the state 
variable of node i by Xi = or Xi — 1, corresponding to 
states S or /, respectively. A pair of nodes share a 
weighted symmetric link Wq > representing their con- 
tact rate. Note that in general these contact rates can 
have any real positive value, unlike network models that 
are based on binary interactions. Susceptible nodes can 
become infected, and infected nodes can then become 
susceptible again upon recovery. This is the well-studied 
SIS (susceptible- infected-susceptible) model [2j. We also 
consider the SIR (susceptible-infected-recovered) model 
where infected individuals become immune to the disease 
upon recovery. Each link is designated either caretaker 
(C) or regular (R), and the fraction of C links is denoted 
p c . We denote this signature of a link by er^ = 1 if the 
link is a caretaker link and <7y = —1 if it is regular. 
These two classes represent different ways of respond- 
ing to an epidemic. Caretaker relationships cause nodes 
to increase their contact frequency Wij if an attached 
node is infected, while regular relationships cause nodes 
to avoid each other (decreasing contact rates). At each 
time step a susceptible i can become infected by one of 
its infected neighbors with a probability that increases 
with link weight. We assume that: 

Pi = 1 - exp (-air) (1) 

where r is the propensity of disease transmission follow- 
ing a contact, and cm = J^. the susceptible's 
contact rate with infecteds. 

An infected individual i recovers with propensity ft 
which yields the probability of recovery 

n = i- cx P (-ft) (2) 

We consider two scenarios: 1) Infected nodes recover at 
a uniform rate ft — j3 or 2) with variable probability. In 
the latter case, caretaker relationships increase a node's 
recovery probability ft according to 



where ft is the base recovery rate, and ft the enhanced 
recovery rate induced by the action of caretakers. The 
quantity Ui represents the total exposure of an infected 
to caretakers and is given by 

*i = ^^Wijil + aij), 

3 
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Fig. 1 (a) An initial network with all nodes susceptible 
(left) has two caretaker links (green) and three regular links 
(black). After the infection of the central node (shown by 
change to red color), regular-linked nodes react by "avoiding" 
the infected node (represented here by increasing distance). 
Caretaker-linked nodes, on the other hand, react by further 
increasing contact rates (represented here by decreasing dis- 
tance), (b) Another network consists of two clusters around 
two central infected nodes (red) . When considering the "care- 
taker effect", the more caretaker interactions (green) a node 
is exposed to, the greater its recovery rate (shown by node 
size; larger nodes have faster recovery rates). Thus after a 
time step, the lower infected node is more likely to recover, 
shown by its transition to susceptible status (blue). 



thus Oi is the total weight of caretaker interactions that 
node i experiences. The parameter ctq sets the scale for 
this exposure. The shape of the sigmoid curve can be 
controlled by the exponent n. 

The infectious state of the system is defined by the 
states Xi of each node. We model the adaptive nature of 
the network weights Wij according to 



StWij = fiUi^Xi + Xj) - 7 (i 



(3) 



Here the first term acts as the driving force of weight 
change, governed by the rate parameter \i. If a link is a 
caretaker link (er^ =1), and one of the adjacent nodes 
is infected (x,i — 1 or Xj — 1), this term is positive and 
causes the weight to increase (if both nodes are infected 
the change is additive). Regular links (ov, = —1), on the 
other hand decrease in strength if one of the connected 
nodes is infected. The second term acts as a restorative 
force, governed by the rate parameter 7 <C fi. Because we 
investigate a system in discrete time we use the following 
update rule for the weights: 

Wij(t+1) = Wij(t) exp [/zo-y (xi + xj) - 7 (wij(t) - w- )] , 

(4) 

a discrete time reformulation of Eq. (|3| . 
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Fig. 2 Infected density (I* = I/N) for SIS dynamics as 
a function of time for different caretaker proportions p c , 
where caretakers do not improve recovery. Erdos-Renyi net- 
works with adaptive rewiring were used (solid lines), as well 
as a similar static network (no rewiring, dashed line). Solid 
lines were obtained by averaging over 100 simulations, so a 
single-simulation plot is overlaid in each adaptive scenario 
for reference. The plot corresponds to Io = 10 2 , N = 10 3 , 
Per = 0.008, fj, = 0.05, 7 = 0.037, j3 = 0.15, r = 0.18. 



3 Results 

We first consider SIS dynamics. At each time step, a ran- 
domly chosen node i can transition from S to I with 
probability pi, or from I to S with probability as 
given above. To study the effect of adaptive rewiring, we 
first consider a system without the caretaker effect on 
the recovery rate, i.e. /?i = (3q. Caretakers only increase 
their interaction with infected individuals. We consider 
a network with weights initially distributed uniformly 
between and 1. Results are shown in Fig. [2] In the 
absence of caretaker links (p c = 0), the equilibrium en- 
demic state I* = It/N is much lower than compared to 
the static network (without rewiring). This is expected, 
as only regular (negative) interactions exist that decrease 
in response to the epidemic. The total network weight 
adapts to a smaller value, decreasing the endemic state. 
The dynamics of the disease and adaptation of the net- 
work is visible in the damped oscillation of the fraction 
of infecteds. 

However, as the fraction of caretakers is increased, 
diseases can attain higher endemic states than their static 
network counterparts. The caretaker dynamics increases 
the interaction rate with infecteds, effectively yielding a 
higher disease prevalence, which is expected. 

The system that lacks a positive caretaker effect rep- 
resents a somewhat artificial limiting case. We therefore 
consider a positive caretaker effect: caretaker relation- 
ships lend higher recovery rates (5\ > /?o to infected 
individuals, see Eq. ^. In particular, we consider the 
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Fig. 3 Extinction probability p ext for SIS dynamics as a 
function of caretaker proportion p c for various values of /3j in 
an Erdos-Renyi network. Note a critical p c value at which ex- 
tinction is maximized. Approaching this value from the left 
yields a gradual increase in extinction, while increasing p c 
past this critical value causes a rapid decrease in p ex t- The 
plot corresponds to Io = 10 2 , N = 10 3 , /i = 0.05, 7 = 0.037, 
r = 0.18, A, = 0.35, a = {a,) I Pba = 0.008. 



effect of varying the maximum recovery rate ft and the 
fraction of caretaker links p c on the extinction proba- 
bility of the disease. The results are depicted in Figs. [3] 
and[4j In general, increasing ft yields higher extinction, 
since caretaker links are more effective at raising recovery 
rates. One would then expect that increasing the care- 
taker proportion p c would also yield higher extinction, 
as more relationships would cause increasing recovery 
rates. However, this is not necessarily the case. Raising 
the caretaker proportion past some ft-dependent critical 
value allows diseases to persist. This critical value also 
serves as a threshold, as increasing p c above this value 
rapidly decreases the extinction probability to 0. This 
is illustrated in Fig. [3j Increasing p c at first yields and 
increased p ext until a maximum is reached. A further 
increase leads to a rapid decrease in extinction proba- 
bility. For the Erdos-Renyi network, the critical fraction 
of caretakers is approximately p c w 10%. For p c values 
above or below this, high extinction probability is seen 
only for very high values of ft. Note however, that even 
for very small fractions of caretakers, a substantial in- 
crease in extinction probability is observed. This suggests 
that, if the caretaker-effect is taken into account, the best 
strategy to extinguish a disease is the existence of a few 
effective caretaker relationships, that safely avoids the 
negative effects that emerge beyond the critical concen- 
tration. Note also that for non-vanishing p c , guaranteed 
extinction (p ex t — 1) is observed only for very high values 
of ft. 

Note that these results were obtained for an Erdos- 
Renyi network. In order to investigate the interaction 



of network adaptation in combination with strong net- 
work heterogeneity, we investigated the dynamics in a 
scale free topology. The results are depicted in Fig. [4] In 
contrast to the Erdos-Renyi system, we observe a high 
extinction of the disease for a wide range of caretaker 
concentrations and recovery parameters ft. The disease 
is endemic in the adaptive, scale free network only for 
small ft and large p c . The implications of these results 
are interesting: In a scale free adaptive network, regular 
links that decrease when connected to infected nodes are 
sufficient to extinguish a disease, even in the presence of 
a considerable fraction of caretaker links. This strongly 
contrasts with the behavior observed in static scale free 
networks, in which the existence of strongly connected 
hubs generally facilitate the spread of a disease. In the 
adaptive network, for p c <C 1, it is sufficient that the 
majority of nodes decrease their interactions with the in- 
fected subpopulation. In scale-free networks, hubs that 
possess a large number of links will adaptively reduce the 
majority of their regular weights, and thus their ability 
to serve as a gateway of the disease to spread throughout 
the network. In this regime, the effect of caretaker rela- 
tionships and their effect on recovery are benign. Only 
when the fraction of caretaker links reaches a large value 
such that also hubs become predominantly caretakers, 
the situation changes, and the disease will evolve into an 
endemic state. 

To explain these results, consider a susceptible node i 
and its total rate of interaction with infected neighbors: 



The ratio of SI interaction rates and total interaction 
rate a = J2i<j w ij is given by 



asi 



= -^5>S/(i)(l-s 4 ) 



Averaging this measure over the time-course of a disease 
gives us a measure of the typical fraction of contacts due 
to SI interaction: 
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Now consider this time averaged (asi) as a function of 
p c for various values of ft, see Fig. p)J For ft = /3 (i.e. 
no caretaker effect on recovery rates), the rate of SI in- 
teractions increase steadily as p c is increased, yielding a 
more stable endemic state and high prevalence. When 
the caretaker effect is taken into account, we observe an 
initial decrease of SI interactions until a critical value 
is reached below which the disease will go extinct, indi- 
cated by the solid line. Increasing p c further can result in 
increasing SI interactions beyond this critical value, en- 
tering a regime in which a large fraction of caretaker links 
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Fig. 4 Two-parameter phase diagrams showing extinction probability for SIS dynamics as a function of maximum caretaker 
effectiveness /3i and caretaker proportion p c . Erdos-Renyi (left) and Scale- Free networks (right) were considered. In the black 
regions, extinction probability is while extinction probability is 1 in the white regions. The plots suggest that increasing 
the caretaker proportion past a critical value yields a decreased extinction probability in both networks. On the Erdos-Renyi 
network, p c w 1CP 1 yields maximum disease extinction, while extinction is most likely for p c « on the Scale-Free network. 



The plots correspond to Io = 10 , N 
Pi = Po, Per = 0.008, (Erdos-Renyi) 



= W 3 , n = 0.05, 7 = 0.037, r = 0.18, Po = 0.35, a = \ t=Q if (a,) \ 
and a mean degree fco = 2 in the scale free network. 



> otherwise 



results in a negative effect. In the scale-free network, the 
qualitative behavior is similar. The crucial difference is 
that typically, the adaptive process of regular links is suf- 
ficient to put the fraction of SI links below the critical 
value even in the absence of caretaker links. 

Next we turn out attention to the effect of caretaker 
adaptive networks on systems that are better described 
by SIR dynamics. Here individuals (nodes) exist in one 
of three states, susceptible (S), infected (I) or recovered 
(R). Individuals can transition from S to I with prob- 
ability pi and from I to R with probability rj, as given 
above in Eqs. ([!]) and The state R is absorbing, so 
once all infected nodes in a population recover, the dis- 
ease dies out (see Fig. [6]). In order to investigate the 
impact of caretaker dynamics and an SIR scenario, we 
focus on the attack rate (ratio) and the epidemic peak. 
The attack ratio (AR) is simply the fraction of the popu- 
lation which contracts the infection at some point during 
the epidemic. Since every infected node eventually enters 
the recovered class, this is equivalent to the fraction of 
recovered nodes at the end of the epidemic: 



AR 



If 



The epidemic peak (EP) is the maximum infected frac- 
tion attained in the population over the course of the 
epidemic. Figure [7] depicts the attack ratio as a function 
of p c for various values of the recovery rate parameter 
P\. Interestingly, without a caretaker effect = (3 ) 
the increase in attack ration is not substantial as p c is 
increased. For fi\ > /3o, we observe a decrease in at- 
tack ratio even for small fractions of caretaker links. The 
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Fig. 6 Infected density (J* = I/N) for SIR dynamics as a 
function of time for different caretaker proportions p c . Erdos- 
Renyi networks with adaptive rewiring were used, as well as 
a similar static network (no rewiring, dashed line). Increasing 
p c lowers the epidemic peak as well as the attack rate. Note 
also that the static network trace closely resembles the p c = 
trace, showing that SIR diseases in this system are not 
significantly affected by dynamic link weights alone. The plots 
correspond to I = 25, N = 10 3 , /x = 0.05, 7 = 0.037, r = 
0.45, Po = 0.20, cr = (ffi) j t=0 if (at) \ t=Q > otherwise 

Pi = Po, n — 2, per = 0.008, (Erdos-Renyi). Scale-Free 
network results were similar. 



minimum attack ratio is attained only in a regime where 
most links are caretaker links. 
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Fig. 5 Time-averaged SI contact fraction (asi) for SIS dynamics with different values of the caretaker proportion p c . Three 
/3i values were chosen, 0.35 (circles), 0.60 (dots), 0.80 (arrows) to correspond with low, intermediate, and high traces in the 
phase diagram of Fig. El An Erdos-Renyi network was used (left), as well as a Scale-Free network (right). The horizontal solid 
lines represent a critical value for (asi) above which the extinction probability vanishes and below which the disease goes 
extinct. The plots correspond to I = 10 2 , N = 10 3 , fi = 0.05, 7 = 0.037, r = 0.18, /3 = 0.35, cr = (ex,) j t=Q if (a t ) \ t=Q > 
otherwise j3i = /So, Per = 0.008, (Erdos-Renyi) and mean degree ko — 2 (Scale-Free). 



Figure [8] depicts the attack ratio as a function of 
both system parameters j3± and p c and compares the be- 
havior in both network architectures, Erdos-Renyi and 
Barabasi-Albert. In contrast with the SIR system, net- 
work topology does not substantially change the dynam- 
ics, both networks exhibit a similar attack ratio as a func- 
tion of /?i and p c . For fixed /?i increasing p c first decreases 
the attack ratio until a minimum is attained. Increasing 
p c further increases the attack ratio again. A consistent 
effect is observed in the response of the epidemic peak 
to changes in fli and p c , see Fig. [9] 

The dynamics seen above for the attack rate are mir- 
rored in the epidemic peak EP as well (Fig. [9]) , which 
decreases as caretaker effectiveness (represented by (3{) 
increases. There is again a critical relationship with p Cl 
as values of p c rs 0.2 tend to minimize the epidemic peak 
for /?i > /3q. Again though, for /3\ — /3q, increasing p c 
yields a monotonic increase in EP. 



4 Conclusions 

Individual response can have a great impact on the dy- 
namics of spreading diseases on complex networks. In 
particular, if one uses an avoidance strategy whereby all 
individuals simply avoid infecteds, the endemic state of 
an SIS disease can be drastically reduced. On the other 
hand, allowing individuals (caretakers) to become closer 
to infecteds is a calculated risk. If the caretakers are 
not effective healers (such as non-physician parents and 
children), then the severity of the disease generally in- 
creases. But if the caretakers are effective healers (con- 
sider doctor/patient relationships, for example), then the 




Fig. 7 Attack rate AR as a function of p c for SIR dynamics 
with various values of /3i in an Erdos-Renyi network. For each 
Pi > Po, the attack rate is minimized for some value of p c 
between 10" 1 and 10°. As /3i increases, this minimum point 
shifts subtly to the right. This shows that the more effective 
caretakers are at healing, the more caretaker relationships 
the system can permit before they have a negative impact on 
the attack rate. The plots correspond to Iq — 25, N = 10 3 , 
H = 0.05, 7 = 0.037, r = 0.25, f3 Q = 0.20, a = (a,) \ t=Q if 
("») | t=0 > otherwise j3 r = /3 , n = 2, p E R = 0.008. 



outcome of the disease is generally improved even by a 
small number of them. If too many caretakers are in- 
troduced, though, their healing benefit is overridden by 
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Fig. 8 Two-parameter phase diagrams showing the dependence of attack rate in SIR dynamics on maximum caretaker 
effectiveness /3i (normalized by the baseline-recovery probability /3o) and caretaker proportion p c . Erdos-Renyi (left) and 
Scale-Free (right) networks were considered. Attack rate approaches zero in the white regions, while it approaches 1 in the 
black regions. Note that increasing p c yields lower attack rates for p c < 0.2, but increasing past this critical value yields 
increasing attack rates. There is a critical value p c « 0.2 at which attack rate is minimized for most values of Pi . Furthermore, 
this effect is seen in both ER and SF networks, though attack rates are lower overall on the SF network. The plots correspond 
to I = 25, TV = 10 3 , fx = 0.05, 7 = 0.037, r = 0.25, /3 = 0.20, a = (<r 4 ) | t=Q if (a t ) \ t=Q > otherwise fa = /9o, n = 2, 
Per = 0.008, (Erdos-Renyi) and mean degree ko = 2 (Scale-Free). 




Fig. 9 Two-parameter phase diagrams showing the dependence of the epidemic peak (EP) in SIR dynamics on the maximum 
caretaker effectiveness /3i (normalized by the baseline-recovery probability /?o) and caretaker proportion p c . Erdos-Renyi (a) 
and Scale-Free (b) networks were considered. The epidemic peak approaches zero in the white regions, while it approaches 1 
in the black regions. Note the similarities to the attack rate diagram in Fig. [8] The epidemic peak is minimized for p c « 0.2 
for most values of but for p c < 0.2 or p c > 0.2, the attack rate is greater for a given value of The plots correspond 
to I = 25, JV = 10 3 , (J, = 0.05, 7 = 0.037, r = 0.40, A) = 0.20, a = (cr 4 ) | t=Q if (a t ) \ t=Q > otherwise fa = /3 , n = 2, 
Per = 0.008, (Erdos-Renyi) and mean degree ko = 2 (Scale-Free). 



their increased exposure, yielding a worse outcome than 
if the population had simply not reacted. 

These findings have a number of implications in pub- 
lic health. For one, in a large-scale epidemic there cer- 
tainly exists a critical fraction of doctors and aid workers 
in the population. If there are too few or too many, they 
can actually increase the total number of individuals in- 



fected over the course of the disease. In such cases, it 
would actually be more beneficial to employ an avoid- 
ance strategy whereby all individuals, including doctors 
and aid workers, simply avoided infected individuals. In 
the particular case of SIS endemic diseases, we have seen 
that the critical caretaker proportion is actually p c = 
on Scale-Free networks. This suggests that networks that 
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exhibit a strong variability in interaction statistics and 
at the same time are adaptive, are less susceptible to 
the risk of endemic diseases, and the natural instinct to 
avoid infection is more effective in eliminating a disease 
than the positive effects that caretakers may have. 
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